A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.).

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

Genetic Mapping, Identification, and Characterization of a Candidate Susceptibility Gene for Powdery Mildew in Cannabis sativa.

Powdery mildew (PM) in Cannabis sativa is most frequently caused by the biotrophic fungus Golovinomyces ambrosiae. Based on previously characterized variation in susceptibility to PM, biparental populations were developed by crossing the most resistant cultivar evaluated, 'FL 58', with a susceptible cultivar, 'TJ's CBD'. F1 progeny were evaluated and displayed a range of susceptibility, and two were self-pollinated to generate two F2 populations. In 2021, the F2 populations (n = 706) were inoculated with PM and surveyed for disease severity. In both F2 populations, 25% of the progeny were resistant, while the remaining 75% showed a range of susceptibility. The F2 populations, as well as selected F1 progeny and the parents, were genotyped with a single-nucleotide polymorphism array, and a consensus genetic map was produced. A major effect quantitative trait locus on C. sativa chromosome 1 (Chr01) and other smaller-effect quantitative trait loci (QTL) on four other chromosomes were identified. The most associated marker on Chr01 was located near CsMLO1, a candidate susceptibility gene. Genomic DNA and cDNA sequencing of CsMLO1 revealed a 6.8-kb insertion in FL 58, relative to TJ's CBD, of which 846 bp are typically spliced into the mRNA transcript encoding a premature stop codon. Molecular marker assays were developed using CsMLO1 sequences to distinguish PM-resistant and PM-susceptible genotypes. These data support the hypothesis that a mutated MLO susceptibility gene confers resistance to PM in C. sativa and provides new genetic resources to develop resistant cultivars. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Functional analysis of Salix purpurea genes support roles for ARR17 and GATA15 as master regulators of sex determination.

The Salicaceae family is of growing interest in the study of dioecy in plants because the sex determination region (SDR) has been shown to be highly dynamic, with differing locations and heterogametic systems between species. Without the ability to transform and regenerate Salix in tissue culture, previous studies investigating the mechanisms regulating sex in the genus Salix have been limited to genome resequencing and differential gene expression, which are mostly descriptive in nature, and functional validation of candidate sex determination genes has not yet been conducted. Here, we used Arabidopsis to functionally characterize a suite of previously identified candidate genes involved in sex determination and sex dimorphism in the bioenergy shrub willow Salix purpurea. Six candidate master regulator genes for sex determination were heterologously expressed in Arabidopsis, followed by floral proteome analysis. In addition, 11 transcription factors with predicted roles in mediating sex dimorphism downstream of the SDR were tested using DAP-Seq in both male and female S. purpurea DNA. The results of this study provide further evidence to support models for the roles of ARR17 and GATA15 as master regulator genes of sex determination in S. purpurea, contributing to a regulatory system that is notably different from that of its sister genus Populus. Evidence was also obtained for the roles of two transcription factors, an AP2/ERF family gene and a homeodomain-like transcription factor, in downstream regulation of sex dimorphism.

Cannabinoids function in defense against chewing herbivores in Cannabis sativa L.

In the decades since the first cannabinoids were identified by scientists, research has focused almost exclusively on the function and capacity of cannabinoids as medicines and intoxicants for humans and other vertebrates. Very little is known about the adaptive value of cannabinoid production, though several hypotheses have been proposed including protection from ultraviolet radiation, pathogens, and herbivores. To test the prediction that genotypes with greater concentrations of cannabinoids will have reduced herbivory, a segregating F2 population of Cannabis sativa was leveraged to conduct lab- and field-based bioassays investigating the function of cannabinoids in mediating interactions with chewing herbivores. In the field, foliar cannabinoid concentration was inversely correlated with chewing herbivore damage. On detached leaves, Trichoplusia ni larvae consumed less leaf area and grew less when feeding on leaves with greater concentrations of cannabinoids. Scanning electron and light microscopy were used to characterize variation in glandular trichome morphology. Cannabinoid-free genotypes had trichomes that appeared collapsed. To isolate cannabinoids from confounding factors, artificial insect diet was amended with cannabinoids in a range of physiologically relevant concentrations. Larvae grew less and had lower rates of survival as cannabinoid concentration increased. These results support the hypothesis that cannabinoids function in defense against chewing herbivores.

Evolution of a ZW sex chromosome system in willows.

Transitions in the heterogamety of sex chromosomes (e.g., XY to ZW or vice versa) fundamentally alter the genetic basis of sex determination, however the details of these changes have been studied in only a few cases. In an XY to ZW transition, the X is likely to give rise to the W because they both carry feminizing genes and the X is expected to harbour less genetic load than the Y. Here, using a new reference genome for Salix exigua, we trace the X, Y, Z, and W sex determination regions during the homologous transition from an XY system to a ZW system in willow (Salix). We show that both the W and the Z arose from the Y chromosome. We find that the new Z chromosome shares multiple homologous putative masculinizing factors with the ancestral Y, whereas the new W lost these masculinizing factors and gained feminizing factors. The origination of both the W and Z from the Y was permitted by an unexpectedly low genetic load on the Y and this indicates that the origins of sex chromosomes during homologous transitions may be more flexible than previously considered.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives.

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

Correlations among morphological and biochemical traits in high-cannabidiol hemp (Cannabis sativa L.).

Cannabis sativa is cultivated for multiple uses including the production of cannabinoids. In developing improved production systems for high-cannabinoid cultivars, scientists and cultivators must consider the optimization of complex and interacting sets of morphological, phenological, and biochemical traits, which have historically been shaped by natural and anthropogenic selection. Determining factors that modulate cannabinoid variation within and among genotypes is fundamental to developing efficient production systems and understanding the ecological significance of cannabinoids. Thirty-two high-cannabinoid hemp cultivars were characterized for traits including flowering date and shoot-tip cannabinoid concentration. Additionally, a set of plant architecture traits, as well as wet, dry, and stripped inflorescence biomass were measured at harvest. One plant per plot was partitioned post-harvest to quantify intra-plant variation in inflorescence biomass production and cannabinoid concentration. Some cultivars showed intra-plant variation in cannabinoid concentration, while many had a consistent concentration regardless of canopy position. There was both intra- and inter-cultivar variation in architecture that correlated with intra-plant distribution of inflorescence biomass, and concentration of cannabinoids sampled from various positions within a plant. These relationships among morphological and biochemical traits will inform future decisions by cultivators, regulators, and plant breeders.

Structural variation of a sex-linked region confers monoecy and implicates GATA15 as a master regulator of sex in Salix purpurea.

The Salicaceae, including Populus and Salix, are dioecious perennials that utilize different sex determination systems. This family provides a useful system to better understand the evolution of dioecy and sex chromosomes. Here, a rare monoecious genotype of Salix purpurea, 94003, was self- and cross-pollinated and progeny sex ratios were used to test hypotheses on possible mechanisms of sex determination. To delimit genomic regions associated with monoecious expression, the 94003 genome sequence was assembled and DNA- and RNA-Seq of progeny inflorescences was performed. Based on alignments of progeny shotgun DNA sequences to the haplotype-resolved monoecious 94003 genome assembly and reference male and female genomes, a 1.15 Mb sex-linked region on Chr15W was confirmed to be absent in monecious plants. Inheritance of this structural variation is responsible for the loss of a male-suppressing function in what would otherwise be genetic females (ZW), resulting in monoecy (ZWH or WWH ), or lethality, if homozygous (WH WH ). We present a refined, two-gene sex determination model for Salix purpurea, mediated by ARR17 and GATA15 that is different from the single-gene ARR17-mediated system in the related genus Populus.

De Novo Assembly and Annotation of 11 Diverse Shrub Willow (Salix) Genomes Reveals Novel Gene Organization in Sex-Linked Regions.

Poplar and willow species in the Salicaceae are dioecious, yet have been shown to use different sex determination systems located on different chromosomes. Willows in the subgenus Vetrix are interesting for comparative studies of sex determination systems, yet genomic resources for these species are still quite limited. Only a few annotated reference genome assemblies are available, despite many species in use in breeding programs. Here we present de novo assemblies and annotations of 11 shrub willow genomes from six species. Copy number variation of candidate sex determination genes within each genome was characterized and revealed remarkable differences in putative master regulator gene duplication and deletion. We also analyzed copy number and expression of candidate genes involved in floral secondary metabolism, and identified substantial variation across genotypes, which can be used for parental selection in breeding programs. Lastly, we report on a genotype that produces only female descendants and identified gene presence/absence variation in the mitochondrial genome that may be responsible for this unusual inheritance.

Composition and functional properties of hemp seed protein isolates from various hemp cultivars.

Hemp seed protein isolates (HPI) were extracted from seven commercial hemp cultivars, a Cornell breeding line, and a commercial hemp heart product, and their composition and functional properties were investigated. HPI contained different ratios of edestin, vicilin, and albumin proteins, which affected protein solubility and functionality. Higher protein solubility was associated with cultivars that contained more vicilin and albumin, which influenced the subsequent functional properties of HPI. Significant differences in water holding capacity (0.83-1.05 g water/g protein isolate), oil holding capacity (1.28-1.81 g oil/g protein isolate), foam capacity (52.9%-84.9%), and foam stability (68.1%-89.4%) were observed across HPI. The Cornell hemp breeding line exhibited the highest protein solubility at pH 7.0 and was uniquely capable of forming an emulsion. The relationship identified between hemp seed protein composition and functionality, in conjunction with the demonstration of an on-going hemp breeding line, suggest that continued, targeted development of hemp cultivars can improve its seed protein functional properties for ingredient utilization in plant-based foods.

Protoplast-Based Transient Expression and Gene Editing in Shrub Willow (Salix purpurea L.).

Shrub willows (Salix section Vetrix) are grown as a bioenergy crop in multiple countries and as ornamentals across the northern hemisphere. To facilitate the breeding and genetic advancement of shrub willow, there is a strong interest in the characterization and functional validation of genes involved in plant growth and biomass production. While protocols for shoot regeneration in tissue culture and production of stably transformed lines have greatly advanced this research in the closely related genus Populus, a lack of efficient methods for regeneration and transformation has stymied similar advancements in willow functional genomics. Moreover, transient expression assays in willow have been limited to callus tissue and hairy root systems. Here we report an efficient method for protoplast isolation from S. purpurea leaf tissue, along with transient overexpression and CRISPR-Cas9 mediated mutations. This is the first such report of transient gene expression in Salix protoplasts as well as the first application of CRISPR technology in this genus. These new capabilities pave the way for future functional genomics studies in this important bioenergy and ornamental crop.

First Report of Downy Mildew caused by Pseudoperonospora cannabina on Cannabis sativa in New York.

Hemp (Cannabis sativa <0.3% tetrahydrocannabinol) is an emerging crop used for grain, fiber, and cannabinoid production (Fike et al. 2020). In New York, hemp is grown both in controlled environment facilities, including greenhouses, and as a field crop. In August 2020, downy mildew-like symptoms were observed on leaves and inflorescence of hemp plants in a field research trial in Ithaca, NY. Several cultivars, including 'Auto CBD', were affected. Disease was severe with some plants reaching 75% disease severity at the individual plant level. In the most severely affected plots, there was no marketable yield. The disease was characterized by chlorotic and necrotic lesions producing sporangiophores under high humidity. Pigmented sporangia were produced on branched sporangiophores. On artificially inoculated leaves incubated at 18°C, 80% humidity, 12h light for 5d, sporangiophores produced 8-19 pigmented, lemon-shaped sporangia with mean ± SD dimensions of 25.2 ± 3.0 (18.9 to 30.4) x 18.2 ± 2.1 (14.6 to 23.2) µm (n=50). Each sporangium produced 2-5 zoospores after less than 45 min in water at room temperature (22°C). Sporangia were collected from sporulating lesions and DNA was extracted as outlined in Crowell et al. (2020). Fragments of the ribosomal internal transcribed spacer (ITS) region (White et al. 1990), the beta-tubulin ras-associated ypt1 gene (Moorman et al. 2002), and the mitochondrial cytochrome B oxidase subunit 2 (cox2) gene (Hudspeth et al. 2000) were amplified by PCR and sequenced bidirectionally. Sequences were deposited in GenBank under accession numbers OK086084, OM867581, and OM867580, respectively. BLAST searches using the amplified ITS and cox2 sequences resulted in 100% identity to Pseudoperonospora cannabina (HM636051.1, HM636003.1) with ypt1 aligning at 97.95% identity (382/390 bp) with P. cannabina (KJ651402.1). The molecular characterization identified the causal agent as P. cannabina. A representative isolate was deposited in the Cornell Plant Pathology Herbarium as CUP-070922. Sporangia were rinsed from detached leaves and used to confirm pathogenicity on whole plants. Ten 4-week-old 'Anka' plants were spray-inoculated until run off with a suspension of 1x104 sporangia mL-1. Ten control plants were sprayed with water. After inoculation, plants were placed in a 19˚C growth chamber with a 12-h photoperiod and misted for 30 min twice daily to maintain humidity above 80%. Sporangia and previously described symptoms were observed 7 days post-inoculation, while control plants were asymptomatic. The pathogen was reisolated onto detached leaves of 'Anka' from inoculated leaves where both sporangia and oospores were observed. The reisolated pathogen was confirmed morphologically and molecularly, through PCR amplification and bidirectional sequencing of the ITS, cox2, and ypt1 genes, as P. cannabina. To our knowledge, this is the first report of P. cannabina causing hemp downy mildew in New York. Depending on the severity and timing of infections, this disease could pose a significant threat to hemp production in the state. Other members of the genus, P. cubensis and P. humuli cause downy mildew on cucurbits and hops, respectively. As these can cause devastating diseases on their hosts, P. cannabina must be monitored with vigilance as an emerging pathogen (Purayannur et al. 2021; Savory et al. 2011). Literature Cited: Crowell, C. R., et al.2020. Plant Dis. 104:2949. DOI 10.1094/PDIS-04-20-0718-RE Fike, J. H., et al. 2020. Page 89 In: Sustainable Agriculture Reviews, vol 42. Springer, Cham, Switzerland. DOI 10.1007/978-3-030-41384-2_3 Hudspeth, D. S. S., et al. 2000. Mycologia 92:674. DOI 10.2307/3761425 Moorman, G. W., et al. 2002. Plant Dis. 86:1227. DOI 10.1094/PDIS.2002.86.11.1227 Purayannur, S., et al. 2021. Mol. Plant Pathol. 22:755. DOI 10.1111/mpp.13063 Savory, E. A., et al. 2011. Mol. Plant Pathol. 12:217. DOI 10.1111/j.1364-3703.2010.00670.x White, T. J., et al. 1990. Page 315 In: PCR Protocols. A Guide to Methods and Applications. Academic Press, San Diego, CA. DOI 10.1016/B978-0-12-372180-8.50042-1.

Identification and mapping of major-effect flowering time loci Autoflower1 and Early1 in Cannabis sativa L.

Flowering time is an important trait for all major market classes of hemp (Cannabis sativa), affecting yields and quality of grain, fiber, and cannabinoids. C. sativa is usually considered a short-day plant, flowering once night length reaches a critical threshold. Variations in flowering time within and across cultivars in outdoor grown populations have been previously identified, likely corresponding to genetic differences in this critical night length. Further, some C. sativa are photoperiod insensitive, colloquially referred to as "autoflowering." This trait has anecdotally been described as a simple recessive trait with major impacts on phenology and yield. In this work, the locus responsible for the "autoflower" trait (Autoflower1), as well as a major-effect flowering time locus, Early1, were mapped using bulked segregant analysis. Breeder-friendly high-throughput molecular marker assays were subsequently developed for both loci. Also detailed are the flowering responses of diverse cultivars grown in continuous light and the result of crossing two photoperiod insensitive cultivars of differing pedigree.

Repeated turnovers keep sex chromosomes young in willows.

BACKGROUND: Salicaceae species have diverse sex determination systems and frequent sex chromosome turnovers. However, compared with poplars, the diversity of sex determination in willows is poorly understood, and little is known about the evolutionary forces driving their turnover. Here, we characterized the sex determination in two Salix species, S. chaenomeloides and S. arbutifolia, which have an XY system on chromosome 7 and 15, respectively. RESULTS: Based on the assemblies of their sex determination regions, we found that the sex determination mechanism of willows may have underlying similarities with poplars, both involving intact and/or partial homologs of a type A cytokinin response regulator (RR) gene. Comparative analyses suggested that at least two sex turnover events have occurred in Salix, one preserving the ancestral pattern of male heterogamety, and the other changing heterogametic sex from XY to ZW, which could be partly explained by the "deleterious mutation load" and "sexually antagonistic selection" theoretical models. We hypothesize that these repeated turnovers keep sex chromosomes of willow species in a perpetually young state, leading to limited degeneration. CONCLUSIONS: Our findings further improve the evolutionary trajectory of sex chromosomes in Salicaceae species, explore the evolutionary forces driving the repeated turnovers of their sex chromosomes, and provide a valuable reference for the study of sex chromosomes in other species.

Evidence of Asexual Overwintering of Melampsora paradoxa and Mapping of Stem Rust Host Resistance in Salix.

Melampsora rust is a devastating disease of shrub willow in North America. Previous work has identified Melampsora paradoxa as one of two identified rust species in New York State that infect Salix purpurea and other important Salix host species, however little is known about the population of this rust species in this region. Genotyping-by-sequencing was used to identify single nucleotide polymorphisms (SNPs) and assess population diversity of M. paradoxa isolates collected from three Salix breeding populations in Geneva, NY between 2015 and 2020. Statistical analyses of SNP revealed that all isolates collected were clonally derived even though they were collected across years. In 2020, isolates were collected from stem infections where uredospore pustules were observed, and these isolates were also identical to M. paradoxa collected in previous seasons. These data suggest that M. paradoxa sampled across multiple years overwintered and reproduced asexually and that stem infection is a possible mechanism for overwintering, both of which are novel findings for this rust species. Additionally, field disease ratings were conducted on a S. purpurea × S. suchowensis F1 breeding population with high disease severity, enabling the discovery of QTL for resistance on chromosomes 1 and 19. Lastly, Colletotrichum salicis was frequently associated with stem rust and may play a role in M. paradoxa stem infection. Together, this work is the first substantial exploration into M. paradoxa population biology, stem infection, and host resistance in Salix.

Genetic mapping of sexually dimorphic volatile and non-volatile floral secondary chemistry of a dioecious willow.

Secondary chemistry often differs between sexes in dioecious plant species, a pattern attributed to its possible role in the evolution and/or maintenance of dioecy. We used GC-MS to measure floral volatiles emitted from, and LC-MS to quantitate non-volatile secondary compounds contained in, female and male Salix purpurea willow catkins from an F2 family. Using the abundance of these chemicals, we then performed quantitative trait locus (QTL) mapping to locate them on the genome, identified biosynthetic candidate genes in the QTL intervals, and examined expression patterns of candidate genes using RNA-seq. Male flowers emitted more total terpenoids than females, but females produced more benzenoids. Male tissue contained greater amounts of phenolic glycosides, but females had more chalcones and flavonoids. A flavonoid pigment and a spermidine derivative were found only in males. Male catkins were almost twice the mass of females. Forty-two QTL were mapped for 25 chemical traits and catkin mass across 16 of the 19 S. purpurea chromosomes. Several candidate genes were identified, including a chalcone isomerase associated with seven compounds. A better understanding of the genetic basis of the sexually dimorphic chemistry of a dioecious species may shed light on how chemically mediated ecological interactions may have helped in the evolution and maintenance of dioecy.

Mapping the sex determination region in the Salix F1 hybrid common parent population confirms a ZW system in six diverse species.

Within the genus Salix, there are approximately 350 species native primarily to the northern hemisphere and adapted to a wide range of habitats. This diversity can be exploited to mine novel alleles conferring variation important for production as a bioenergy crop, but also to identify evolutionarily important genes, such as those involved in sex determination. To leverage this diversity, we created a mapping population by crossing 6 Salix species (Salix viminalis, Salix suchowensis, Salix integra, Salix koriyanagi, Salix udensis, and Salix alberti) to common male and female Salix purpurea parents. Each family was genotyped via genotyping-by-sequencing and assessed for kinship and population structure as well as the construction of 16 backcross linkage maps to be used as a genetic resource for breeding and selection. Analyses of population structure resolved both the parents and F1 progeny to their respective phylogenetic section and indicated that the S. alberti parent was misidentified and was most likely S.suchowensis. Sex determining regions were identified on Salix chromosome 15 in the female-informative maps for seven of the eight families indicating that these species share a common female heterogametic ZW sex system. The eighth family, S. integra × S. purpurea, was entirely female and had a truncated chromosome 15. Beyond sex determination, the Salix F1 hybrid common parent population (Salix F1 HCP) introduced here will be useful in characterizing genetic factors underlying complex traits, aid in marker-assisted selection, and support genome assemblies for this promising bioenergy crop.

Comparative transcriptomics and eQTL mapping of response to Melampsora americana in selected Salix purpurea F2 progeny.

BACKGROUND: Melampsora spp. rusts are the greatest pathogen threat to shrub willow (Salix spp.) bioenergy crops. Genetic resistance is key to limit the effects of these foliar diseases on host response and biomass yield, however, the genetic basis of host resistance has not been characterized. The addition of new genomic resources for Salix provides greater power to investigate the interaction between S. purpurea and M. americana, species commonly found in the Northeast US. Here, we utilize 3' RNA-seq to investigate host-pathogen interactions following controlled inoculations of M. americana on resistant and susceptible F2 S. purpurea genotypes identified in a recent QTL mapping study. Differential gene expression, network analysis, and eQTL mapping were used to contrast the response to inoculation and to identify associated candidate genes. RESULTS: Controlled inoculation in a replicated greenhouse study identified 19 and 105 differentially expressed genes between resistant and susceptible genotypes at 42 and 66 HPI, respectively. Defense response gene networks were activated in both resistant and susceptible genotypes and enriched for many of the same defense response genes, yet the hub genes of these common response modules showed greater mean expression among the resistant plants. Further, eight and six eQTL hotspots were identified at 42 and 66 HPI, respectively. The combined results of three analyses highlight 124 candidate genes in the host for further analysis while analysis of pathogen RNA showed differential expression of 22 genes, two of which are candidate pathogen effectors. CONCLUSIONS: We identified two differentially expressed M. americana transcripts and 124 S. purpurea genes that are good candidates for future studies to confirm their role in conferring resistance.

Nonadditive gene expression is correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow (Salix spp.).

Many studies have highlighted the complex and diverse basis for heterosis in inbred crops. Despite the lack of a consensus model, it is vital that we turn our attention to understanding heterosis in undomesticated, heterozygous, and polyploid species, such as willow (Salix spp.). Shrub willow is a dedicated energy crop bred to be fast-growing and high yielding on marginal land without competing with food crops. A trend in willow breeding is the consistent pattern of heterosis in triploids produced from crosses between diploid and tetraploid species. Here, we test whether differentially expressed genes are associated with heterosis in triploid families derived from diploid Salix purpurea, diploid Salix viminalis, and tetraploid Salix miyabeana parents. Three biological replicates of shoot tips from all family progeny and parents were collected after 12 weeks in the greenhouse and RNA extracted for RNA-Seq analysis. This study provides evidence that nonadditive patterns of gene expression are correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow. Expression-level dominance was most correlated with heterosis for biomass yield traits and was highly enriched for processes involved in starch and sucrose metabolism. In addition, there was a global dosage effect of parent alleles in triploid hybrids, with expression proportional to copy number variation. Importantly, differentially expressed genes between family parents were most predictive of heterosis for both field and greenhouse collected traits. Altogether, these data will be used to progress models of heterosis to complement the growing genomic resources available for the improvement of heterozygous perennial bioenergy crops.

The Melampsora americana Population on Salix purpurea in the Great Lakes Region Is Highly Diverse with a Contributory Influence of Clonality.

Shrub willows (Salix spp.) are emerging as a viable lignocellulosic, second-generation bioenergy crop with many growth characteristics favorable for marginal lands in New York State and surrounding areas. Willow rust, caused by members of the genus Melampsora, is the most limiting disease of shrub willow in this region and remains extremely understudied. In this study, genetic diversity, genetic structure, and pathogen clonality were examined in Melampsora americana over two growing seasons via genotyping-by-sequencing to identify single-nucleotide polymorphism markers. In conjunction with this project, a reference genome of rust isolate R15-033-03 was generated to aid in variant discovery. Sampling between years allowed regional and site-specific investigation into population dynamics, in the context of both wild and cultivated hosts within high-density plantings. This work revealed that this pathogen is largely panmictic over the sampled areas, with few sites showing moderate genetic differentiation. These data support the hypothesis of sexual recombination between growing seasons because no genotype persisted across the two years of sampling. Additionally, clonality was determined as a driver of pathogen populations within cultivated fields and single shrubs; however, there is also evidence of high genetic diversity of rust isolates in all settings. This work provides a framework for M. americana population structure in the Great Lakes region, providing crucial information that can aid in future resistance breeding efforts.